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, The energy release in gamma-ray bursts is one of the most useful clues on the 

\ nature of their "inner engines". We show that, within the framework of the relativistic 

"^H \ external shocks afterglow model, the narrowness of the observed X-ray luminosity of 

Qh| GRB afterglows, implies that the energy of the GRB jets after the early afterglow phase 

O ' spans less than one order of magnitude. This result is not affected by uncertainties in 

' the electrons energy, magnetic field strength, and external medium density. We argue 

■ that the afterglow kinetic energy is within a factor of two of the initial energy in the 

J> . relativistic ejecta, therefore the energy output of the central engine of long duration 

^ , GRB has an universal value. 

' Subject headings: gamma-rays: bursts - ISM: jets and outflows - radiation mechanisms: non- 

thermal - shock waves 

In the last three years we have learned a great deal about long duration gamma-ray bursts 
(GRBs). The Italian-Dutch satellite BeppoSAX provided angular position of several dozen long 
bursts to within about 3 arc-minutes which enabled follow up observations in the x-ray (see Piro, 
2000), optical, milli-meter and radio frequencies which has provided a wealth of information on 
these explosions. These observations are described well by the relativistic fireball model (see e.g. 
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Piran, 1999). According to this model the energy from the central source is deposited in material 
that moves with speed very close to the speed of hght. The kinetic energy of this material is 
converted to the observed electromagnetic radiation as a result of collisions between fast moving 
material that catches up with slower moving ejecta, and the shock heated circum-burst medium. 
The nature of the "inner engines" that expels relativistic material, which is responsible for the the 
GRBs, is not determined yet. The energy of the relativistic matter ejected by the "inner engine", 
Erei, is one of the most important clues on its nature. Our goal is to find a reliable estimate of 

Erel- 

Given an observed 7-ray fiuence and the redshift to a burst one can easily estimate the energy 
emitted in 7-rays, E^^iso assuming that the emission is isotropic. E^^iso can also be estimated from 
the BATSE catalogue by fitting the flux distribution to theoretical models (Cohen & Piran, 1995; 
Schmidt, 2001). As afterglow observations proceeded, alarmingly large values (Kulkarni et al. 1999) 
(3.4x lO^^ergs for GRB990123) were measured for E^^iso- However, it turned out (Rhoads, 1999; Sari 
Piran k. Halpern, 1999) that GRBs are likely beamed and E^^iso would not then be a good estimate 
for the total energy emitted in 7-rays. We define instead: E^ = [O"^ /2)E^^iso- Here 9 is the effective 
angle of 7-ray emission, which can be estimated from ^5, the time of the break that appears latter 
in the afterglow light curve (Rhoads, 1999; Sari Piran & Halpern, 1999): 9 = {).l2{n/ E^^iYl^tl'^^y^, 
where £'51 is the isotropic-equivalent energy kinetic energy during the adiabatic fireball phase, 
discussed below, in units of lO^^ergs, and ti^days is the break time in days. Recently Frail et al. 
(2001, hereafter Fl) estimated Ej for 18 bursts, finding typical values around lO^^ergs. While 
Ej is closer to Erei it is still not a good estimate. First, we have to take an unknown conversion 
efficiency of energy to 7-rays into consideration: Erai = e~^E^ = e'^{9'^ /2)Ej^iso. Second, the large 
Lorentz factor during the 7-ray emission phase, makes the observed E^y rather sensitive to angular 
inhomogeneities of the relativistic ejecta (Kumar & Piran, 2000). 

We consider here another quantity: Ex.adi the kinetic energy of the ejecta during the adiabatic 
afterglow phasc^. Clearly: E^^i > E^ + Ex, ad = ^-^rei + E^^ad, where E'^, is the angular average 
of E~^. The inequality arises from possible energy losses during the early afterglow radiative phase. 
However observations of long time tails of GRBs suggest that, unless this energy is radiated at 
extremely high energy channel, this losses arc not large (Burenin et al, 1999, Giblin et al, 1999; 
Tkachenko et al., 2000) . Therefore, with e ~ 10% (Kobayashi, Piran & Sari, 1997) (e cannot be 
too close to unity otherwise there won't be afterglow) we expect that Ex^ad ~ -E're/ to within better 
than a factor of 2. Hereafter we drop the subscript kin denoting E = Ex^ad- 

The purpose of this paper is to determine the spread of E using the x-ray afterglow flux. The 
advantage of the method presented here is that it is independent of the uncertain density of the 
ISM, and in fact all other parameters, except for the observed width of the distribution of the jet 
opening angle. 



'The external relativistic shock becomes adiabatic about 1/2 an hour after the explosion, and furthermore the 
loss of energy during the earlier radiative phase is typically not large. 
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One way of determining how E is distributed is to figure out the energy for individual bursts 
by modeling their afterglow emission over a wide range of frequency and time. This procedure, 
carried out for 8 well studied bursts (Panaitescu & Kumar 2001, hereafter PKOl), gives the mean 
energy to be ~ 3 x lO^'^ergs and the standard deviation of the log of energy distribution to be about 
0.3. However, the detailed modeling of individual GRB afterglow emission is cumbersome and time 
consuming, and hard to carry out for a large sample of bursts especially considering that we need 
data in radio, optical and x-ray bands with good time coverage for this kind of an analysis to be 
useful. 

Moreover, this procedure is not necessary if we simply want to know the width of energy distri- 
bution. For this we can use the x-ray afterglow flux at a fixed time after the explosion. The width 
of the distribution of this flux, an easily measurable quantity, yields the width of the distribution 
for the energy release in the explosion. This method is described below. The observation of the 
x-ray flux should be carried out at a sufficiently late time such that the angular variations and fluc- 
tuations across the surface of the ejecta have been smoothed out. This occurs several hours after 
the explosion when the bulk Lorentz factor of the ejecta has decreased to about 10 at which time 
we see a good fraction of the rclativistic ejecta. Conveniently, this is also when the observed x-ray 
in the 2-10 keV band is above the cooling frequency in which case the observed flux is independent 
of the density of the medium in the vicinity of the explosion (Kumar 2000). Furthermore, it is 
best to carry out observations while F > 9~^, that is before the jet begins its sideways expansion, 
which makes the interpretation of the observed flux much simpler. In five cases of GRB afterglow 
light curves where we see the effect of the finite opening angle of explosion as steeping of the light 
curve, we find the effect manifests itself at least one day after the explosion. Thus, the above two 
requirements suggest it is best to consider the x-ray afterglow flux between several hours and a day 
after the explosion. 

The x-ray afterglow fluxes from GRBs have a power law dependence on v and on the observed 
time t (Piro, 2000): fu{t) oc v^^t^" with a ~ 1.4 and /3 ~ 0.9. The observed x-ray flux per unit 
frequency, fx, is related, therefore, to, L^, the isotropic luminosity of the source at redshift, z by: 

W = (TtI) + '^^'"^ ^ Ut)z{z) , (1) 

where Z{z) is a weakly varying function of z. For bursts with 0.5 < z < 4 and with /3 — a ~ —0.5 
we find az ~ 0.31 (for a cosmology with = 0.3 and fl\ = 0.7). Here and thereafter we denote 
by ax the standard deviation of the log(X), unless noted otherwise. 

Assuming that the x-ray luminosity does not evolve with redshift we can relate the dispersion 
of log{Lx) at a fixed observer time after the explosion, aL^, with the observed dispersion af^: 
al^ = aj^ +al ^ aj^ . Using 21 BeppoSAX bursts (Piro, 2000) we find a/^ ^ 0.43 ± 0.1 (see 
the caption for Figs. 1 and 2 for the details of the observations and the analysis), and therefore, 
c"L^ ~ 0.43 to within 25%. This result is supported by 10 x-ray light-curves of GRBs with known 
red-shifts, and a and (3. 
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The x-ray flux, in an energy band above the cooHng frequency at a fixed time after the burst, 
depends on the energy per unit soHd angle in the explosion (provided that T > at the time 
of observation), and on the fractional energy taken up by electrons, eg. The flux does not depend 
on the density of the surrounding medium, n, or its stratification or the fractional energy in the 
magnetic field, es (Kumar, 2000). The flux has a weak dependence on n when the electron cooling 
is dominated by the inverse Compton scattering (Panaitescu &; Kumar 2000). 

Under these, rather general, conditions the standard synchrotron fireball model implies that 
the isotropic equivalent fiux at frequency u above the cooling frequency, at a fixed elapsed time 
since the explosion, is given by (Kumar 2000; Preedman &; Waxman, 2001): 

(p+2)/4 



Lx = Vp 



dE 

In 



-1 (p-2)/4 



(2) 



where dE/dCl is the energy per unit solid angle, and rjp is a constant". Assuming that there is 
no correlation between the microscopic variables, eg, e^, p and dE/dO, we obtain from the above 
equation that CFdE/dCi < (^L^- Using gl^ ~ 0.43 it 0.1 for the 21 BeppoSAX bursts we find that 
(^dE/dQ < 0.43 ±0.1. 

Prom (T^E/dQ can now obtain aE that characterizes the distribution of the kinetic energy 
provided we 

know ao using the trivial relation: cr^E/dn ~ ^E '^^e- Panaitescu & Kumar (PKOl) and 
Frail ct al. (FOl) have estimated the jet opening angles from the observed (or lack of) breaks in 
the light curves of optical afterglow light curves. For 8 GRBs from the PKOl sample we have: 
ag « 0.31 ± 0.06 while a sample of 10 bursts from FOl yields ag « 0.28 ± 0.05. If these values 
are representative for the whole GRB population we find a marginally viable solution within two a 
errors of aE < 0.2 (for the PKOl result) and aE < 0.27 (for the FOl data); to get a viable solution 
we had to take both the values of aL^ one standard deviation above the mean and the value of 
ag one SD below the mean. This result suggests that there is a narrow energy distribution; the 
FWHM of E being less than a factor of 5. If E and 9 are correlated the above relation is modified 
i.e. c'^E/dSi ~ '^E ^'^e ~ ^C<^E<^e- Both the PKOl and the FOl data show that this correlation (Q 
if non-zero is weak, less than 0.35. With such a correlation the allowed energy distribution could 
be somewhat (but not significantly) broader. 

A stronger constrain on aE can be obtained using = 0.3 for 8 GRBs analyzed by PKOl. 
It follows from Eq. 2 that a non-zero value for a^^ makes the distribution for dE/dO, and hence E 
even narrower, however to quantify this effect we need a larger data set. 

We have argued before that E = Ex,adi discussed here, is a rather good estimate to Ej-ei the 
total energy emitted by the "inner engine". The constancy of E^ad is another indication for it 



"it is worth emphasizing that equation (2) is independent of the details of the jet structure i.e., the variation of 
Lorentz factor across the jet, since at about 1/2 a day after the explosion the jet Lorentz factor has dropped to ~ 10 
and we see a good fraction of the entire jet surface. 
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being a good measure of Erei- The constancy of Ejc^ad is also an indication that the assumptions 
that have lead to Eq. 2 are justified. Otherwise it would have been remarkable if starting from 
different levels of initial energy and having different amounts of energy losses the final kinetic energy 
of the afterglow would converge to a constant value. At present there is no way to tell whether 
significant amounts of additional energy is released in other forms, e.g. non relativistic particles or 
neutrinos. Similar arguments suggest further that this total energy emitted by the "central engine" 
does not vary significantly and that it is rather close to the energy estimated here. 

The distribution of relativistic kinetic energy during the afterglow phase is narrow, with full 
width at half maximum less than one decade. These results suggest that the wide distribution of 
directly and indirectly determined E^^iso results from the distribution of beaming angles, from a 
variation of dE/dQ across the jet, and from a variable efficiency in conversion of kinetic energy 
to 7-rays. The fact that GRB engines are "standard" engines in terms of their energy output 
provide a very severe constraint on the nature of these enigmatic explosions. For instance, in the 
collapsar model for GRBs the central engine is composed of a black hole (BH) and an accretion 
disk around it (Woosley 1993; Paczynski, 1998; MacFadyen &; Woosley, 1999). This model has two 
energy reservoirs which can be tapped to launch a relativistic jet: the BH rotation energy and the 
gravitaional energy of the disk. Our result of nearly constant energy in GRBs implies that the 
mass accretion on to the BH plus the possible conversion of rotational energy of the BH to kinetic 
energy of the jet does not vary much from one burst to another inspite of the fact that both the 
disk mass and the BH spin are expected to vary widely in the collapse of massive stars. 

BeppoSAX is a mission of the Italian Space Agency (ASI) with participation of the Dutch 
space agency (NIVR). TP acknowledges support by the US-Israel BSF. 
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Fig. 1. — The distribution of X-ray fluxes (2-10 koV) at t=ll hours after the GRB in 21 afterglows observed by 
BeppoSAX. The sample includes all the fast observations performed by BeppoSAX on GRB from January 1997 to 
October 1999. Data are from Piro (2001, proc. of X-ray astronomy 99, and references therein) Stratta et al. (2001, 
in preparation), Nicastro et al. 1999 (IAUC7213), Feroci et al. (2001 A&A, in press), Montanari et al. (2001, Proc. 
of the 2nd wokshop on GRB in the afterglow era, Rome, 17-20 Oct. 200, E. Costa, F. Frontera J. Hjort eds., in 
press), Piro et al. (1999, GCN 409) and f Zand et al. (2000, ApJ, 545, 266). No X-ray afterglow was detected in 
GB990217 (Piro et al. 1999, lAUC 7111) to the hmiting instrumental sensitivity of 10~" erg cm"^ s~\ 6 hours 
after the burst. In the case of GB970111 (Feroci et al. 1998 A&A 332, L29) a candidate was detected, but evidence 
of fading behaviour is marginal, so we have considered both cases as upper limits (indicated by the arrow in fig) . 
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Fig. 2. — Likelihood contour lines (corresponding to 99%, 90% and 69% confidence levels) in the log(/a;), 
a/^ plane for the X-ray flux distribution as inferred from 21 GRBs detected by BeppoSAX from Jan- 
uary 1997 to October 1999. We determine log(/a;) and af^ by minimizing the likelihood function S = 
— In I [27r((T| -|- cr|^)] exp [— (log — log fx)^ /2{af + a^^)]^; where log /a;; and ai are the observed x-ray 
flux 11 hr after the onset of the i-th GRB and the associated measurement error respectively. The maximal likelihood 
is at log(/^) = -12.2±0.2 and tj/^ = 0A3t {l. Two upper hmits of 10~" orgs cm"^ s"^ in the 2-10 kcv band at 11 
hours after the bursts are included in this data set. The value of cr/^ is 0.42 if we exclude these upper limits. We 
have checked a posteriori with a Kolgomorov-Smirnov test that the distribution is consistent with a gaussian (at 90% 
confidence level). We also note that the predicted number of X-ray afterglows with a flux lower than about 2 x 10~^^ 
is 3.5, consistent with the observed number of objects. 



